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X-RAY  TOPOGRAPHIC  MEASUREMENT  OF  STRAIN  FIELDS 


S.R.  Stock*,  Haydn  Chen  and  H.K.  Birnbaura 

Department  of  Metallurgy  and  Mining  Engineering 
University  of  Illinois  at  Urbana-Champaign 
Urbana,  IL  61801  USA 


INTRODUCTION 

X-ray  diffraction  topography  offers  a  unique  combination  of 
advantages  for  imaging  dislocations  and  accumulated  plastic 
deformation  in  single  crystals.  Its  strain  sensitivity  and 
penetrating  power  are  greater  than  that  of  electron  microscopy,  and 
topography  can  be  used  non-destructively  to  examine  a  single 
specimen  many  times  during  the  course  of  an  experiment.  These 
attributes  combine  with  the  limit  in  attainable  magnlf ication  to 
make  topography  most  useful  for  characterizing  crystals  with  low 
dislocation  densities.  With  the  advent  of  synchrotron  radiation 
sources  and  of  rapid  imaging  systems  for  laboratory  sources,  the 
emphasis  is  shifting  from  characterization  studies  to  dynamic,  in- 
situ  observation  of  experiments.  In  heavily  deformed  crystals, 
equi-incllnation  contours  (analogues  to  bend  contours  in  electron 
microscopy)  have  been  noted'  '  although  no  quantitative  use  of  the 
angular  information  implicit  In  a  set  of  contours  has  been  made. 

The  purpose  of  this  paper  is  to  use  sets  of  equi-inclination 
contours  (EIC)  to  determine  the  components  of  strain  around  a  strain 
center  about  which  individual  defects  cannot  be  resolved.  Analytic 
methods  are  developed  for  both  the  monochromatic  and  synchrotron 
white  radiation  sources.  Some  preliminary  results  are  reported  for 
a  precipitate  of  p-NbH  in  a  niobium  single  crystal. 


*Now  at  the  Dept,  of  Mat.  Sci.  &  Eng.,  Northwestern  Univ. ,  Evanston, 
IL  61201,  USA. 


Reproduced  from 
be«t  eveileble  copy. 


RELATIONSHIP  BETWEEN  THE  STRAIN  TENSOR  AND  THE  ROTATION  REQUIRED  FOR 
DIFFRACTION  AT  DIFFERENT  POINTS  IN  A  CRYSTAL 


Consider  a  strain  center  which  deforms  a  thin,  rtbbon-like 
single  crystal.  A  rectilinear  coordinate  system  with  Its  origin  at 
the  strain  center  may  be  defined  with  normal  to  the  largest  face 
of  the  crystal  (Fig.  1).  Diffraction  from  planes  normal  to  axis  x*, 
denoted  by  P^,  will  be  considered  for  the  case  of  a  parallel, 
spatially-broad  monochromatic  beam  of  x-rays.  Displacement  of  atoms 
from  their  positions  prior  to  the  deformation  Is  given  by  u^  which 
are  related  to  the  strain  tensor  by  e,  ,  =  l/2(5u,/ax,  +  hu./ax.). 
Variation  in  the  interplanar  spacing  or1Jin  the  orientation  Jof  the 
diffracting  planes,  will  limit  diffraction  to  a  portion  of  the 
specimen.  Only  deformation  with  components  normal  to  P^  will  be 
sensible,  however,  and  only  changes  in  the  displacement  will 
change  the  diffraction  conditions.  The  change  in  the  displacement 
across  a  unit  volume  element,  du^/Ax  is  related  to k  the  angular 
change  in  diffraction  conditions,  J  i.e.  to  Ato^i.A+l) ,  the 
rotation  required  to  align  a  new  region  for  diffraction.  The 
superscript  k  defines  the  rotation  axis  x^»  the  subscript  1  gives 
the  diffraction  plane  P^  and  the  arguments  l  and  A+l  identify  the 
spatial  positions  of  the  two  volumes  of  material  aligned  for 
diffraction.* 


<  a  )  (  b  > 


Fig.  1  Schematic  diagram  of  the  equl-inclinatlon  contour 
analysis,  (a).  The  contour  position  on  a  topograph  is 
shown  at  distance  X.  from  the  strain  center  for  a  rotation 
of  from  the  undeforraed  region  of  the  single  crystal. 

(b).  The  resulting  relationship  between  X^  and  w  is  plotted 
for  results  obtained  from  topographs  taken  with  diffraction 
vectors  h^  and  2h^. 


Since  the  derivation  of  the  equation  relating  Au^  and  du^/SXj 


*The  rotation  axis  x^  is  that  normally  used  in  x-ray  topography;  xi 
and  2L3  diffraction  from  P3  and  P^,  respectively,  in  transmission 
and  either  xj^  or  353  for  diffraction  from  P2  in  reflection. 


is  the  same  for  diffraction  from  each  P^,  only  that  from  Pj  will  be 


explicitly  considered.  Of 

dilational  strain  and  bu^/hx 
strains.  The  total  rotat 


u.  the  terms  bu^/ax, ,  bu^ /ax.  ... 
j/dx„  and  ^Uj/bx^  are  components  or  s 
otatlon  Aojj,  is,  therefore,  composed 
separate  contributions  from  the  variation  in  e^ch  of  the  terms 
du./9x,.  The  form  of  the  contributions  to  Acu^  are  well  known * 
but  the  presentation  of  the  particular 
necessary  to  understand  the  analysis. 


is  a 
shear 
of 


2) 


formulism  used  here  is 


The  change  in  the  Bragg  angle  A9..  from  the  varying  dilational 
field  is  related  to  the  change  in  bu^/ox^  from  position  £  to  £+1 
through  the  differential  form  of  Bragg's  law: 


Ad/d  =  -A9_cot9 
O  Do 


Aei;l  =  AdUj/ax^Jl.Jl+l), 


(1) 


where  Ad  =  (d  j-dQ)  “  (d  -dQ)  and  dQ,  d  and  d^+1  are  the  d- 
spacings  far  "trora  the  strain  center  (i.e.  "in  the  undeformed  state), 
at  £  and  at  £+1,  respectively;  A9D=  0£+j-0.  anc*  ®o’  and 
the  Bragg  angles  of  the  undeformed  crystal,  at  £  and  at  £+1, 
respectively. 


£+1 


are 


The  shear  terms,  Qu^/ax^  and  au^/bx^  represent  tilting  about 
axes  2L3  and  respectively.  With  T:he  experimental  rotation  axis 
jc-j,  tne  relationship  of  the  terms  and  au^/ax^  to  the 
rotation  required  to  orient  region  £+1  for  diffraction  will  be  quite 
different.  The  contribution  of  au^/bx-  is  trivial  to  calculate: 
the  angle  A9g^  through  which  the  crystal  must  be  rotated  is  simply 
the  difference  in  the  amount  of  tilting  at  two  positions: 

A0S1(£,£+1)  (rad.)  =  b^/bx^i+l)  -  aUl/ax2(£)  =  Aau1/ax2.  (2) 

Limitation  of  space  precludes  detailed  consideration  of  the  geometry 
required  to  relate  the  variation  in  bu^/hx^  to  the 
rotation  A9_„;  the  details  will  appear  elsewhere.  Here  we  will  only 
note  that  the  complication  arises  because  the  axis  about  which 
changes  in  bUj/bu-j  occur  is  x>>  while  the  axis  of  rotation  used  to 
produce  EIC  in  the  experiment  is  x^.  The  equation 

sin  9q  =  sin  (0q+  A9g2)  cos  y  (3) 

relates  A0g2  a°d  y=  Aauj/dx^(£,£+l) .  Before  writing  the  expression 
for  At0j(£,£+1),  it  is  informative  to  consider  the  magnitude  of  the 
rotations  arising  from  each  terra.  Assuming  that  each  of  the  deriva¬ 
tives  a^/bXj  has  a  value  of  0.01  and  that  0Q  ■  12. A®  (^  =  [200]  for 
niobium  with  Mo  radiation),  the  rotations  A0y,  A0gl  and  A9g2  are 
450,  2100  and  2.2  seconds  of  arc,  respectively.  The  contribution 
of  Abuj/ax^  is  so  small  that  it  may  be  ignored,  and  one  may  write 

Au^  (rad)  =  -Abu^ /bx^ ( £ , £+1)  tan  0Q  +  AbUj/bx2(£,£+t)  . 


« 


L 


(A) 


This  expression  nay  be  generalized  for  diffraction  from  and 

rotation  about  rotation  axis  jx^,  leading  to 

(rad)  =  -Adu^/dx^i.l+l)  tan  0q  +  Abu^bx^  (A,Jf+l)  (5) 

with  the  subscript  j  denoting  the  third  reference  direction  jr.  which 
is  perpendicular  to  both  2^  and  x^. 

X-RAY  TOPOGRAPHY  AND  THE  MEASUREMENT  OF  THE  COMPONENTS  OF 
STRAIN 

In  x-ray  topographs  of  elastically  deformed  crystals  or  those 
containing  an  excess  dislocation  density  of  one  sign,  diffraction  of 
monochromatic  x-rays  produces  narrow  bands  of  darkening,  termed 
equi-inclination  contours.  By  multiple  exposure  of  a  single  piece 
of  film  or  by  a  set  of  singly  exposed  emulsions,  the  position  of  the 
equi-lncllnatlon  contour  may  be  mapped  as  a  function  of  rotation. 
Characteristic  radiation  or  monochromatized  synchrotron  radiation 
may  be  used  to  produce  the  equi-inclination  contours.  Similar 
contours  (termed  absorption  edge  contours)  may  be  produced  with 
synchrotron  white  radiation  by  orienting  the  deformed  crystal  so 
that  the  range  of  wavelengths  selected  encompasses  that  of  the 
absorption  edge  of  an  element  of  the  specimen  crystal. 

The  relationship  between  equi-inclination  or  absorption  edge 
contours  to  strains  (Eq.  5)  is  discussed  using  Fig.  1  which 
schematically  shows  a  number  of  contours  recorded  between  the  strain 
center  and  the  portion  of  the  crystal  to  which  the  strain  field  does 
not  extend.  Since  transmission  topography  averages  the  strain  over 
the  thickness  of  the  specimen,  only  the  two-dimensional  variation 
in  strain  is  measured,  i.e. ,  bu^/dx .(x^ ,x^) •  Given  constant 
rotation  between  contours,  it  is  simple  to  determine  the  cumulative 
rotation  £rom  the  undeformed  region  to  the  position  of  the  £th 
contour  u^(Ji).  The  separation  between  each  contour  and  the  strain 
center  X.  is  measured  along  _x^,  the  direction  normal  to  the 
diffraction  planes  P^  if  the  usual  rotation  axis  x^  is  used.*  In 
the  example  shown  in  ^Fig.  la  for  diffraction  from  Pj,  and  rotation 
about  the  axis  u)^(X)  ai^d  X.  may  be  measured  for  each  contour 

along  each  line  x=*  c  and  w,  may  be  plotted  as  a  function  of 
(Fig.  lb).  By  taking  the  differences  with  respect  to  the  position 

of  zero  strain,  Eq.  5  becomes 

“i(WC>  “  -3u1/3x1(x1  ,x^-c)  tan  0^+  bui/bx1(xi ,xk»c) .  (6) 

*One  must  allow  for  foreshortening  of  the  diffracted  image  which  can 
be  significant.  One  does  not  need  to  correct  for  the  differences  in 
diffraction  angle  at  different  contours,  however,  as  the  shift  in 
contour  positions  from  this  source  is  insignificant. 


r 


The  individual  terms  bu^/bx.  and  bu^/bx.  cannot  be  separated  based 
solely  upon  Eq.  6  and  a  single  plot  of  cumulative  rotations  vs. 
position. 

The  required  additional  information  may  be  obtained  by  carrying 
out  the  above  procedures  for  the  first  and  higher  order  reflections 
from  diffraction  planes  (represented  by  diffraction  vectors  hj 
and  2hj).  The  two  rotation  vs.  distance  curves  which  result  due  to 
the  difference  in  tan  0  are  shown  as  solid  and  dashed  lines  in  Fig. 
lb  for  measurements  taken  with  h^  and  2h^,  respectively.  The 
difference  in  cumulative  rotation  at  between  the  two  measurements 
corresponds  to 

k  k  hj  2hj 

«l<xifXk-c>|  -  wi(xi  ’3VC)  I  =  au1/ex1(K1>xk=c)(tan0o  -tan0Q  1)(7) 
2h^  h^ 

allowing  determination  of  bu^/bx  (x.pX^c)  and  bu^/bx. (x.,xk=c) . 
Values  of  bu  /bx^x^.x^)  and  bu^bx  (x^x,  )  may  be  folind1  by  repeated 
use  of  Eq.  7  for  contours  intersecting  various  x^.  This  approach 
allows  the  following  terms  to  be  determined  as  a  function  of 
position  (x|,x^):  for  diffraction  plane  Pj  and  rotation  axis 
bu./bx^  and  bUj/bx2  are  measured;  for  diffraction  plane  and 

rotation  axis  x i »  bu^/bx^  and  bu^/bx2  are  measured  and  for 

diffraction  plane  P2  and  rotation  axis  Xj  or  jc-j,  bu2/bx2  anc* 
bu^/dx^  or  bu2 /bx^  are  measured. 

The  strain  tensor  is  not  completely  defined  by  these  measure¬ 
ments,  however,  as  bu^bx.  and  bu^/bx.  remain  to  be  determined.  For 
a  thin-  specimen,  plane  stress  conditions  apply  but  add  no  informa¬ 
tion  which  helps  determine  bu^/bx  and  bu^/bx  .  Use  of  additional 
rotation  axes  are  required  if  these  terras  are*  to  be  measured.  For 
diffraction  from  Pj  (required  for  determination  of  bu  /bx^),  the 
appropriate  rotation  axis  is  Xg.  Terra^  bUj/bXj  and  bUj/bx2  contri-. 
bute  negligibly  to  the  rotation  Awj ,  and  1  only  bUj/bx^  has  a 
significant  effect.  For  bu^/bx^,  diffraction  from  P^  and  rotation 
about  Xg  *8  required;  bu^/bx^and  ou2/bx2  produce  no  contribution. 

Instead  of  using  diffraction  from  first  and  higher  order 
planes,  one  may  use  two  or  more  wavelengths  and  a  single  diffraction 
vector  for  the  analysis.  This  method  is  particularly  useful  when 
tunable  synchrotron  radiation  is  utilized.  Another  alternative  is 
to  use  diffraction  vectors  h^  and  -h.,  i.e.  exchange  entrance  and 
exit  sides  of  the  specimen.  The  dilat tonal  and  shear  components 
will  produce  rotations  which  will  add  in  one  case  and  be  of  the 
opposite  sense  in  the  other.  Measurements  of  cumulative  rotation 
vs.  distance  for  different  entrance  sides  (diffraction  vectors  h^ 
and  -hf)  result  in  a  difference  given  by 

.  •  h , 

Aw. (h. )  -  Aw. (-h, )  -  -2A  bu,/bx,tan  0  . 

11  11  1  l  o 


(8) 


Fig.  2  A  composite  of  topographs  tak.ec  with  j*  =  [002]  and 
with  the  crystal  rotated  140”  about  axis  [110]  between 
exposures.  Plates  0  and  OM  show  a  multiple  exposure  topo¬ 
graph  of  the  equi-inclination  contours  and  an  optical  micro¬ 
graph  of  the  hydride  positions  in  the  crystal,  respectively. 


Fig.  3  Cumulative  rotation  as  a  Fig.  4  Values  of  some  of  the 

function  of  distance  from  the  0-NbH  bu^/Ox.  along  [110]  from  the 

precipitate  center.  The  contours  centerJof  the  macroscopic 

were  obtained  for  diffraction  from  p-NbH  precipitate. 

[110]  and  rotation  about  [001]. 


MEASUREMENT  OF  STRAIN  TENSOR  COMPONENTS  ABOUT  A  p-NIOBIUM  HYDRIDE 
PRECIPITATE 


The  Nb-H  system  provides  a  paradigm  for  hydrogen  related 
failure  via  stress  assisted  nucleation  and  fracture  of  hydride 
precipitates.  Precipitation  of  g-NbH  and  the  accommodation  of  the 
volume  increase  has  been  investigated  by  electron 
and  with  optical  and  scanning  electron 
The  contour  analysis  described  above  is  suited  to 
determination  of  the  deformation  near  the  hydrides  and  will  be 
applied  below  to  obtain  some  of  the  components  of  strain. 

Hydride  precipitates  were  grown  during  cooling  from  room 
temperature  in  a  thermal  gradient  in  a  NbHQ#Q3  single  crystal  having 
a  thickness  of  75  pm.  Many  small  hydride  precipitates  formed  at  the 
colder  end  while  larger,  well  spaced  precipitates  grew  in  the  warmer 
portions  of  the  crystal.  The  hydride  precipitates  remained  on 
warming  to  room  temperature  because  of  the  considerable  hysteresis 
in  the  reversion'  .  Lang  topographs,  taken  with  Mo  K  radiation 
were  used  to  select  an  isolated  precipitate  for  the  analysis. 

Equi-inclination  contours  about  the  precipitate  were  obtained 
using  reflection  and  transmission  topography  for  rotations 

about  [110]  and  [001]  axes.  The  equi-inclination  contours  in  the 
£  =  [002]  topograph  of  Fig.  2,  plate  0,  were  produced  by  multiple 
exposure  with  rotations  about  the  [110]  axis  of  210  seconds  of  arc 
between  contours.  Deflection  of  contours  between  neighboring 

precipitates  clearly  shows  the  long  range  interaction  of  the  strain 
fields  which  appear  to  extend  two  or  three  millimeters  from  the 
hydride  plates.  The  appearance  of  interacting  contours  is  somewhat 
similar  to  photoelastic  contours  observed  about  model  precipitates 
with  a  similar  volume  expansion'5'.  It  is  notable  that  precipitates 
with  a  large  number  of  contours,  particularly  B  and  D,  are  those 
which  have  a  pair  of  parallel  dendritic  arms  and,  hence,  greater 
volume  expansion.  Figure  2,  _a  composite  of  single  exposure 

topographs  with  rotation  axis  [110]  and  diffraction  vector  - 
[002],  shows  the  positions  of  the  contours  at  140  arc  second 
increments  of  rotation.  The  magnitude  of  deformation  at  precipitate 
B  is  clearly  evident;  the  area  beside  the  precipitate  is  diffracting 
at  angles  greater  than  0.3°  from  the  orientation  at  which 
diffraction  would  occur  if  no  hydride  were  present.  The  sense  of 
rotation  for  the  deformation  field  of  the  precipitate  is  opposite  to 
that  of  the  macroscopic  curvature  of  the  crystal. 

Figure  3  shows  the  cumulative  rotation  vs.  distance  from  the 
center  of  the  precipitate  for  equi-inclination  contours  about  the 
midpoint  of  precipitate  B  for  the  rotation  axis  [001].  The  rotation 
angles  are  measured  relative  to  the  angle  at  which  the  contours  from 
the  precipitate  merge  with  contours  reflecting  the  general  macro¬ 
scopic  curvature  of  the  crystal.  This  position  is  the  point  at 


accompanying 
microscopy^3 
microscopy' 3 


« 


J 


which  the  contributions  from  the  precipitate  and  the  uniform 
curvature  of  the  entire  crystal  cancel.  The  position  corresponding 
to  zero  strain  from  the  precipitate  is  that  at  which  the  uniform 
bending  contours  are  no  longer  deflected  by  the  field  of  the 
precipitate.  Extrapolation  to  zero  strain  is  necessary  to  obtain 
the  magnitude  of  the  strain  component  du^/dx.  at  a  given  position. 
An  estimate  of  the  strain  between  the  last  contour  and  the  reference 

) 

Fig.  4 

presents  estimates  of  some  of  Su^/3x.  along  the  line  _x^  =  0.  The 
results  of  Fig.  4  were  based  on  contour  positions  on  opposite_sides 
of  precipitate  B  taken  with  the  same  diffraction  vector,  =  [110]; 
and  the  shear  and  dilational  terms  contribute  as  if  the  measurements 
were  for  ±  h^,  on  the  same  side  of  the  precipitate. 


position  is  obtained 
(measured  between  the 


by  extrapolating  ( Adu^/dx. )/(A3u. /3x 
e  last  two  contours)  to  ^zero  strai-Ji. 


A  major  difficulty  in  the  analysis  near  hydrides  results  from 
the  complex  hydride  shapes;  it  is  difficult  to  identify  the  "center" 
of  the  deformation.  While  the  indistinct  double  image  of 
precipitate  B  in  the  £  =  [002]  topographs  corresponds  closely  with 
the  pair  of  parallel  dendritic  arms  seen  in  optical  and  scanning 
electron  micrographs,  it  is  very  difficult  to  identify  the  gross 
precipitate  morphology.  Therefore  the  origin  from  which  to  compare 
the  equi-inclinatlon  contours  from  different  reflections  is 
difficult  to  determine.  The  use  of  white  radiation  and  absorption 
edge  contours  will  minimize  this  difficulty  since  images  are 
obtained  from  the  entire  crystal.  Fig.  5  shows  a  set  of  absorption 
edge  contours  of  a  lightly  deformed  niobium  crystal  which  was 
produced  with  diffraction  vector  j[_  =  [002]  and  90  arc  seconds 


change,  from  a  slightly  deformed  niobium  crystal.  (Data 
were  taken  at  Daresbury  Synchrotron  Radiation  Source  in 
collaboration  with  S.T.  Davies  and  D.K.  Bowen.) 
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rotation  between  topographs^\  The  exposure  of  the  Ilford  L4 
nuclear  emulsion  plates  was  chosen  so  that  either  side  of  the 
absorption  edge,  heavily  or  lightly  absorbing,  could  be  reproduced 
in  good  contrast. 

The  inter-contour  spacings  from  any  single  reflection  are  much 
more  reliable  than  comparisons  between  reflections  since  the 
reference  position,  that  of  the  precipitate  image,  will  vary  little 
for  the  rotations  used.  This  can  be  seen  from  Fig.  3  where  there  is 
little  deviation  of  the  individual  points  from  the  average  curves. 
Because  the  difference  of  the  displacement  distance  derivatives 
determine  the  contour  spacings,  the  uncertainty  of  the  origin  from 
reflection  to  reflection  makes  the  absolute  magnitudes  bu^/dx,  some¬ 
what  less  reliable. 

As  described,  the  contour  measurements  may  be  used  to 
characterize  elastic  strains.  For  plastically  deformed  specimens, 
the  technique  is  only  sensitive  to  the  excess  dislocation  density, 
not  to  the  total  dislocation  density. 
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